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Abstract: In the framework of a generalized iterative scheme introduced previously to account for 
the non-analytic coupling dependence associated with the renormalization-group invariant mass scale 
A, we establish the self-consistency equations of the extended Feynman rules (A-modified vertices of 
zeroth perturbative order) for the three-gluon vertex, the two ghost vertices, and the two vertices of 
massless quarks. Calculations are performed to one-loop-order, in Landau gauge, and at the lowest 
approximation level (r = 1) of interest for QCD. We discuss the phenomenon of compensating poles 
inherent in these equations, by which the formalism automatically cancels unphysical poles on internal 
lines, and the role of composite-operator information in the form of equation-of-motion condensate 
conditions. The observed near decoupling of the four-gluon conditions permits a solution to the 2- 
and-3-point conditions within an effective one-parameter freedom. There exists a parameter range in 
which one solution has all vertex coefficients real, as required for a physical solution, and a narrower 
range in which the transverse-gluon and massless-quark propagators both exhibit complex-conjugate 
pole pairs. 




1 Summary of the extended iterative scheme 


The spontaneous emergence of a renormalization-group (RG) invariant mass scale A ]!]] from the renor¬ 
malization process is arguably the most important nonperturbative effect in strictly renormalizable 
quantum field theories, since this quantity sets the scale for all dimensionful observables (except those 
dominated by heavy extraneous masses). Its coupling dependence, 
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(where R denotes a renormalization scheme and u the arbitrary renormalization scale within R, and 
where > 0 in an asymptotically free theory) is non-analytic in a way that will always remain 
invisible in a perturbation expansion around g 2 = 0. Moreover, the several known obstructions || to 
the existence and uniqueness of a Borel transform of the perturbation series all have to do, at least 
qualitatively, with the presence of this scale. There exists, therefore, the intriguing possibility that 
accounting systematically for the A dependence of correlation functions may be the minimal step beyond 
perturbation theory needed to define a strictly renormalizable theory uniquely. (Here ”A dependence” 
does not, of course, refer to a mere reparametrization of the perturbation series, as obtained by 
solving O for g in terms of A, or equivalently by leading-logarithms resummation. Genuinely 
nonperturbative A dependence, exemplified by the way vacuum condensates occur in operator-product 
expansions (OPE), is typically polynomial or inverse-polynomial). 

The present paper elaborates on a specific scheme, outlined earlier in this journal Q, of accounting 
systematically for the A dependence, under the restriction (which in an asymptotically free theory turns 
out to be a weak one) that the known standard technique of renormalization remain applicable with 
at most inessential modifications. This scheme takes the form of an extended iterative solution to the 
integral equations for correlation or vertex functions, starting from a set of extended Feynman rules 
for the superficially divergent basic vertices of the theory, as distinct from the ordinary Feynman rules 
(bare vertices) rd°' pert , whose iteration generates the perturbative series. These extended vertices 
are quantities of zeroth order (p = 0) with respect to the perturbative g 2 dependence but contain a 
“seed dependence” on the nonanalytic scale (1.1), approximated systematically. Under the combined 
requirements of globality (in order to be applicable in loop integrals, the approximation must in 
principle be valid over the entire momentum range) and of the preservation of power counting (as 
a basic prerequisite of standard renormalization technique), the choice of approximating functions 
is remarkably unique: they must be functions rational with respect to A and therefore (since A 
is dimensionful) also with respect to momentum variables. For the representation of the (p = 0) 
nonperturbative A dependence, rational approximants perform the same basic role that polynomial 
approximants play for the perturbative g 2 dependence. The poles of these rational functions, together 
with the numerator zeroes, will in general form discrete approximations to branch cuts of the vertices 
in their complex-momentum planes. We follow ||] in denoting the extended rules by T^’ 0 !, with r the 
(denominator) degree of rational approximation. At each level r, the r[ r,0 l approach the Y^)v ert in 
the “perturbative limit” A —> 0; on the other hand their sequence as r increases may be viewed as an 
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analytic continuation-through-resummation of the zeroth-order terms of the OPE. 

The central problem of such a method is to demonstrate dynamical self-consistency of these ex¬ 
tended Feynman rules. The integral equations (Dyson-Schwinger, or DS, equations) for the set T of 
proper vertex functions Tat, where N stands summarily for the number and types of external legs, 
take the schematic form 

= T { ° )pert + (£-) 2 T7v[r (0)pert ,r], ( 1 . 2 ) 

where go is the bare gauge coupling, and Tat a set of nonlinear dressing functionals, containing loop 
integrals over combinations of T’s. With the inhomogeneous terms always given by the bare p(°)P er ' t ) 
iterating instead around a nonperturbatively modified set of starting functions rl r ’°] can work only 
if the interaction terms 3 >at, at each order l of the iteration, not only generate an Z-th order, ’’quasi- 
perturbative” power correction in g 2 (v), but also reproduce the nonperturbative parts, r[ r,0 l — r(°) peri , 
of the new zeroth-order input: as compared to perturbation theory, the solution must be able to 
’’establish its own zeroth order”. To have the functionals <J>Ar, in spite of their g q prefactor, produce 
certain terms of zeroth order in g 2 is not trivial. Moreover, for the method to be practical, the number 
of extended Feynman rules should remain finite (and small), as in perturbation theory, and in view of 
the notorious infinite hierarchical coupling in eqs. d) - with each 4 >at coupling to T’s up to Tat+i or 
even Tat +2 this is again nontrivial. 

The mechanism discussed in || for simultaneously ensuring both objectives exploits the structure 
(0) of the scale A in conjunction with the renormalizable divergence structure of the theory. To 
briefly describe its main line, let 

{ c }^ = { c r,l) c r,2) • • • c r,fc r } (1-3) 

• • • r qi 

be the complete set of k r numerator coefficients of the level-r rational approximants Ty for all 
superficially divergent vertices (of which, we recall, there are seven in covariantly quantized QCD), 
and let 

{d} [r] = {d r ,l, dr,2, ■ ■ ■ dr } m r } (1.4) 

be the complete set of m r denominator zeroes (pole positions) in units of A 2 . Both c’s and d's are 
dimensionless, real numbers. Then 

r M =r (;)P«-< + A M({ c }M i { d }M;A), (1,5) 


in a notation suppressing momenta and all other variables not immediately pertinent to the argument. 
Here {c|^ denotes the subset of {c}M appearing in the vertex Tat, etc. Upon evaluating, say, the 
first iteration (one-loop order, l = 1) of eq. (L2), dimensionally regularized in D = 4 — 2e, with the 
functions (ED as input, one obtains after some algebraic decomposition, 
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Here {CjM is a set of nonlinear algebraic expressions in the input coefficients ( |1.3| / |1.4|) , while {dy^,^ N 
denotes the subset of denominator roots ( |1.4| ) in the vertices Tjv' other than T^ to which provides 
coupling. The appearance of the first term on the r.h.s. of ( |1.(^ ) is nontrivial: it occurs only if all seven 
basic vertices are treated by mutually consistent, nonperturbative approximants of the same level r. 
This term appears with a prefactor, 
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which, by virtue of an exact RG identity, is independent of the renormalized coupling g{y), and finite 
as e —► 0: 


n ( e ) = ^-[l + 0(e,elne)], 
Po 


independent of g 2 


( 1 . 8 ) 


Here Po is the leading beta-function coefficient of eq. O- Note how in this exact result one coupling 
factor (< 7 q) and one divergence factor Q) get ’’eaten” to produce a coupling-independent and finite 
quantity. 

It therefore becomes possible to reproduce analytically the nonperturbative part, Ay, of the 


zeroth-order input (1.5) by imposing the matching or self-consistency conditions, 
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which says that all basic vertices must exhibit one common set of denominator zeroes (still differing, 
however, for different types of external legs, or basic fields), and 


4cf l ({c}M,{d}M)=<v 


(z — 1,... Ay), 


( 1 . 10 ) 


which ensures reproduction of numerator structures. 

It is crucial that the nonperturbative terms establish themselves in a finite manner, since in this 
way one avoids the introduction of nonlocal counterterms, and thus preserves another basic element 
of standard renormalization. It is equally crucial that the above mechanism, as shown by the last two 
factors of (1.7), is tied to the loop divergences of the integral equations: this gives the superficially 
divergent vertices a privileged position, such that formation of nonperturbative Atv’s remains rigor¬ 
ously restricted to these vertices. In spite of the infinite hierarchical coupling, the number of extended 
Feynman rules does not proliferate, and in fact remains the same as for the bare vertices. 

In what follows, we focus exclusively on this self-consistency process for the generalized Feynman 
rules, and therefore refer the reader to [Sf] for what needs to be said about the last term of eq. 


(1.6) - representing the p = 1 quasi-perturbative correction - and about the perturbative “boundary 
condition” and essentially standard renormalization procedure it requires. We only note, for later use, 
that the condition of having the remaining divergence exactly equal to the perturbative one would 
require to be of the form 
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with z$ the one-loop coefficient of (g/Air) 2 /e in the perturbative renormalization constant Zn for the 
Vertex T jy - a condition which may impose extra constraints on {c}M that for low r may be satishable 
only within approximation errors. 

The exploration of this extended-iterative scheme represents a calculational program of some 
length. It was begun in pj] with an illustrative derivation of eq. (0) and of matching conditions 
o at r = 1 and l = 1 for the transverse two-gluon vertex. While calculations of the ghost and 
fermion two-point functions follow essentially the same pattern, those for the higher superficially di¬ 
vergent vertices, with N = 3 and 4, are not straightforward extensions to more kinematical variables. 
As a result of the subtle interplay between the coupled DS equations, they reveal a whole array of 
new aspects and intricacies. We therefore plan to present and discuss the r = 1, l = 1 self-consistency 
calculations for these vertices in several parts. In the present paper, we focus on the remaining super¬ 
ficially divergent vertices of the gauge, ghost, and massless-quark sectors up to N = 3: the three-gluon 
vertex T3y, ghost vertices T GG and r G y G , where V and G label vector (=gluon) and ghost external 
legs, and fermion vertices T FF and T^y^.. In the companion paper [Q] we will deal with the four-gluon 
vertex T/iy, the highest superficially divergent vertex, which is particularly complicated both kinemat¬ 
ically and in its DS equation. In these two papers, consideration of the fermion (quark) functions will 
be restricted - as were the N = 2V calculations of [3] - to the case of massless quarks, where the 
fermionic mass scales, too, are simply multiples of A. For massive fermions, the presence of ’’extra¬ 
neous” RG-invariant mass scales not having the structure ( |1.1| ) causes additional complications with 
which we plan to deal separately. 

In section 2 of the present article, we recall the DS equation for the T 3 y vertex. While summarizing 
known material, this section seems necessary to establish notation and a precise starting point for the 
subsequent discussion. In the present program we deal exclusively with the ’’ordinary” DS equations, 
without additional Bethe-Salpeter resummations in their interaction terms, in which the distinguished 
leftmost external line always runs into a bare vertex. Section 2.2 focuses on the phenomenon of 
compensating poles in the 3V equation where they make their first appearance. While at first sight a 
merely technical point, these turn out to be an important structural element, by which the formalism 
automatically prevents the appearance of ’’wrong” poles on internal lines. A systematic account of 
these leads to the rearranged integral equation of section 2.3, whose terms now exhibit an extended- 
irreducibility property. It is the rearranged equations that form the most convenient framework for 
the self-consistency problem of the extended Feynman-rules. Extraction of the matching conditions 
( |1.9| - |l-li|) at the lowest level of rational approximation of interest for QCD (r = 1) and at one loop 
(l = 1) is discussed in sect. 2.4 for the 3-gluon vertex. Section 3 considers the equations for the two 
ghost vertices, and section 4 the equations for the two fermion vertices in the massless case. 

The combined 2-plus-3-point self-consistency system is discussed in sect. 5 . A noteworthy result, 
which will continue to hold after inclusion of 4-gluon-conditions, is that the set of denominator coeffi¬ 
cients (0), while restricted by (0), is not fully determined by the divergent parts of DS loops, and 
that composite-operator information, in the form of equation-of-motion condensate conditions ( DS 
equations at coincident spacetime points ) is required at this point to complement the usual equations. 
On the other hand, the system is found to nearly decouple from the 4-gluon one, so that a solution 
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without 4-gluon equations is possible with an effective one-parameter freedom. The results, when 
compared to the more restricted and heuristic attempt of refs. [|[] for a pure-gluon theory, will be 
seen to represent significant progress, and in particular to include the existence of a parameter range 
where one of the several solutions to the nonlinear system has all zeroth-order vertices entirely real in 
the Euclidean, as required for a physical solution. 

As in m, all calculations are performed for the Euclidean theory, and in Landau gauge. The 
Landau gauge provides some welcome reduction of the considerable complexity of loop computation 
with the extended Feynman rules: here the two ghost vertices turn out to remain perturbative, and 
calculations can be restricted to amplitudes with only transverse (if any) gluon legs, which then form 
a closed DS problem. This has the obvious disadvantage that nothing can be inferred as yet about 
the approximate saturation or violation of Slavnov-Taylor identities, which we recall are statements 
about amplitudes with at least one longitudinal gluon leg, and which in the low orders of an iterative 
scheme are not, of course, expected to be exactly self-consistent. However, since the physical degrees 
of freedom of the gauge field are in the transverse sector, one may expect the essential parts of the 
nonperturbative structure to develop here (as is obviously true for the 2-gluon function). Also, we 
do not yet use the most general color-and-Lorentz-tensor structure of the 41/ vertex, which would 
lead to calculations of prohibitive length, but restrict our study to a theoretically motivated tensor 
subset capable of dynamical self-consistency. Nor do we consider the quasi-perturbative corrections, 
(tf) ^ ^iv ^ ' rl notation of eq. s While all these questions are interesting in themselves, they 
must form subjects of future study. 

To the extent that the present scheme uses the DS equations as a framework, its purpose is not to 
provide exact numerical solutions of DS equations at low levels of decoupling; for a review of the work 
in this direction the reader is referred to [jS]. Here the aim is to develop an analytic approximation 
method that provides some insight into the nonperturbative coupling structure, and in particular to 


identify the precise mechanism, connected with the divergence structure, by which the scale (1.1) 
establishes itself in correlation functions. In particular, such a scheme allows qualitative changes in 
the elementary propagators - the appearance of zeroth-order, finite, real or complex mass shifts - to 
be followed in a more transparent fashion. 


2 The three-gluon vertex equation 

2.1 Integral equation and input 

The DS equation for the proper three-vector vertex 5^3 V in Euclidean momentum space is written 
diagrammatically in Fig. 1. The form shown is a compact but hybrid one: most terms on the r.h.s. 
have not been resolved down to the level of proper vertices, but feature connected and amputated 
functions T' which are one-particle irreducible (1PI) only in the horizontal channel of the diagram, 
while otherwise still containing reducible (1PR) terms. Thus in term (A) 3 of Fig. 1, the four-gluon T 
matrix T^ v for the horizontal channel ( denoted there as T’ s ) is to be decomposed further as in Fig. 2: 

Tav = Tav - A\ ( 2 . 1 ) 
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Tky — Ai + A 2 + A 3 + r 4y. 


(2.2) 


Here T 4 y is the proper, fully 1PI, four-gluon vertex, while A\ and A 2 , A 3 are dressed one-gluon re¬ 
ducible terms in the horizontal channel and the two crossed channels, respectively. Analogous relations 
apply to the ghost-antighost-gluon-gluon and quark-antiquark-gluon-gluon T' matrices, Tqgvv anc ^ 
T' F p vv , of terms (B )3 and (E )3 respectively. 

The ’’standard” form of the F 3 y equation in Fig. 1 displays the characteristic asymmetry, common 
to all DS equations, of having the leftmost external leg always ending in a bare vertex, while the other 
legs run into dressed vertices. This structure is at the core of a problem plaguing all treatments (and 
not just the present approximation method) of vertices with N > 3: while the exact solution of the 
equation may be known to have a certain Bose or Fermi symmetry, the equation does not display this 
symmetry manifestly, and approximate solutions to it therefore usually fail to exhibit the full desired 
symmetry. Enforcing the symmetry by imposing extra conditions on the vertex coefficients leads to 
overdetermination in the self-consistency equations. In the framework of an iterative solution, ’’trivial” 
symmetrizations could of course be used to cure this problem, but at the expense of depriving oneself 
of an important indicator of the overall error at the level of approximation considered. 

Since we will be working throughout at the one-loop (l = 1) level, characterized by a single ID- 
dimensional momentum integration in the dressing functional d> 3 y, the term (D )3 of Fig. 1 with two 
DS loops does not yet contribute (its contribution to the l.h.s. of eq. ( 1.10 ) will be of order (1 //?o) 2 )- 

The input for the self-consistency calculation must consist of the Euclidean extended Feynman 

r ip qi 

rules ry J , at the same level r of rational approximation, for all seven superficially divergent vertices. 
With the exception of rj^, whose r = 1 form will be detailed in ||], these have been listed in Q. 
Here we need to recall only two elements carrying special restrictions. First, the gluon-propagator 
rule, 1 #’°! = -(r^ 0] ) _ 1 , will be simplified throughout by adopting the Landau (£ = 0) gauge fixing. 
Then 


{D lrfi] W)% = 5 ab t^(k)Dp°\k 2 ), 
Wk u 


t^(k) = - 


k 2 


where, at the r = 1 level, 


D F ’°\k 2 ) = k 2 + u\ t 1 A 2 + 


Ul,3 


A 4 


-1 


k 2 + ui^A 2 . 
k 2 + ui )2 A 2 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 


(■ k 2 + 0-1,1 A 2 ){k 2 + cry 3 A 2 ) * 

Second, the general color structure of the T 3 y vertex itself, whose self-reproduction we examine in 
this section, 


(r 3 y) a6c — ifabc^'(f) T d abc E(d)i 

will be simplified from the outset to a pure f abc structure, i.e. one puts 


(2.7) 


r (d) ~ 0 


( 2 . 8 ) 
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and omits the / on r tfy The reason is that in the much more complicated color structure of T±y 
discussed in j|], we will disregard those color-basis tensors that would feed the d a bc portion through 
the (A) 3 term of Fig. 1. It is conceivable that some d a b c structure could be made to self-reproduce 
through the other one-loop terms of Fig. 1 alone, but since we view the seven basic DS equations as 
an interrelated whole, it does not seem consistent to us to keep one source of such terms and neglect 
the other. The tensor structure then is 


(T 3V (pi,P2,P3)) P a lc = l fabc{ ^ a (P2-P3) P F 0 (pl,p 2 3 -,p 2 1 ) 

+ 5 ap (p 3 -pi) K F 0 (pl,p 2 1 ;pl) 

+ 5 pK (p 1 - P 2 YFq^pI-pI) 

+ {P2-P3) p {P3-PiY{pi -P2) ct F 1 (p 2 1 ,p 2 2i pI) 

+ [ 10 terms not contributing to totally transverse vertex ] >. (2-9) 


At level r = 1 (and only at r = 1), the invariant functions Fq 1 ' 0 " 1 and F 1 ] 1,0 ^ are conveniently written in 
a form fully decomposed into partial fractions, 


f^ 0 ] (pIp 2 2;pI) =i + zi ) 1 (n 1 + n 2 ) + (s ll 2 + ^)n 1 n 2 + xi ) 3 n 3 
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(aq, 4 + ^) n i + ( X M + 


x l,4 
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n 3 + xi ) 5 nin 2 n 3 , 


( 2 . 10 ) 
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(p 2 i,p 2 2,pI) 
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A 2 


^1,6 (niii2 + n 2 n 3 + n 3 iii) + xi^niirir 


featuring the building blocks 


( 2 . 11 ) 
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A 2 

Pi + u 1,2 A2 


(* = 1,2,3). 


( 2 . 12 ) 


Here we anticipate that the DS self-consistency conditions ( |1.9| ) will enforce u\ 2 = ^ 1 , 2 , i-e., a common 
denominator factor in all gluonic vertices. The factorized (with respect to the 3 variables p 2 ) denomi¬ 
nator structure of these approximants may be viewed as arising from a triple-spectral representation, 


Tb,i = 




dz\dz 2 dz 3 


P0,l(zi,Z2,Z 3 ,pj,pl,pf) 


, 2 \ ’ 


(zi - p()(z 2 - pi)(z 3 - pi) 


(2.13) 


(where the spectral functions are still allowed some polynomial p 2 dependence) by a discrete approxi¬ 
mation. At r = 1, it exhibits a single denominator zero, p 2 = — «y 2 A 2 , in all three variables. 

Although the approximation of branch-cut structures by poles is an old technique, a casual look 
might still suggest a danger here that these poles in vertices could somehow take on a life of their own 
and roam the formalism as unphysical particles. This does not happen, for two reasons. First, DS 
self-consistency will transfer these poles down to the two-point vertices (negative-inverse propagators), 
so that the propagators themselves will develop zeroes at these pole positions, as seen in (^1]). The 
unamputated, connected Green functions, which are the quantities having physical interpretation as 
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propagation amplitudes, will therefore be nonsingular at these positions in the squared momenta of 
their external lines, and no S-matrix elements for unphysical particles of masses ui^A 2 will arise, as 
emphasized already in [:3j]. Second, and perhaps more remarkably, the formalism will also automatically 
cancel the poles of type ( [ 2 . 12 ) when they arise on certain internal lines, so that there will be no 
Cutkosky discontinuities corresponding to production of such objects. To demonstrate this in the 
simplest context will be the subject of the next subsection. 


2.2 ” Compensating” poles in Y iV 


The phenomenon of automatic cancellation of superfluous poles on internal lines seems to have been 
noted first by Jackiw and Johnson and by Cornwall and Norton []7j in 1973. While there are obvious 
technical differences between their (Abelian) models and the present QCD study ( in particular, their 
poles arise in longitudinal-vector channels, whereas here they appear in the transverse-gluon sector ), 
this mechanism is very much the same in both cases. 

Consider the T 3 y equation of Fig. 1 in order [r, 0], and compare residues of both sides at the poles 
in the variable p 2 of the leftmost leg. In a partial-fraction decomposition with respect to p 2 , the l.h.s. 
can be written 


r M] 

1 3V 


i) 

n =1 


A 2 

Pi + Ur,2n^ 2 


+ bIi\p 3 , Pi ) 


+ [ terms with (pi) 1 , (p 2 2 ) 2 ,... (p|) r ]. 


(2.14) 


The compact notation suppresses all tensor structure. All terms displayed still have rational structure 
in p|, p\. The terms in the second line are allowed Q as long as only the conditions of overall 
asymptotic freedom and preservation of perturbative power counting are imposed, but will in fact 
turn out to be more strongly restricted. The residue at p 2 = —u ri 2 nA 2 for the l.h.s. is then A 2 B n 
(n = 1,... r). 

For the r.h.s., to keep the argument as simple as possible, we invoke for the moment all simplifica¬ 
tions available for our specific calculation according to subsect. 2.1, i.e. we disregard terms (D) 3 , (E) 3 , 
and (B) 3 (the general case will be outlined in subsect. 2.3). Then the only denominator structure with 
respect to p\ must come from the T'^y* amplitude of term (A) 3 , for which p 2 = (p 3 +P 1) 2 = se is the 
(Euclidean) Mandelstam variable in the horizontal channel. According to general structure theorems 
on correlation functions ||, residues at poles in such a variable must factorize with respect to the two 
sides of the channel, or be at most sums of factorizing terms. Thus if q±, q 3 are the momenta of the 
loop gluons of term (A) 3 , the T^ v amplitude in the vicinity of se = — 2 nA 2 must behave as 


rp! M] 
1 4V 


1 i / n(Ql,Q3)^n(p3,Pl) 
SE + «r,2nA 2 


+ [ regular terms ]. 


(2.15) 


The partial Bose symmetry of T' iV implies that both residue factors must be given by the same function 
or column vector of functions, ’F ri . Note that ( |2.15| ) in no way contradicts the 1PI property of T^ v 
in the horizontal channel: the pole factor is not a propagator of any of the elementary fields of the 
theory, so the term is technically allowed to appear ( as would, e.g., a bound-state pole ) not only 









iii T' iV but in fact in the fully 1PI piece, of eq. (T2). However, the observation does suggest a 
natural enlargement of the notion of reducibility, as will be discussed below. 

The residue comparison shows that must be proportional to B n , 


^n(P3,Pl) = M n B n (p 3 ,pi), 


(2.16) 


with M n some matrix, and that (again in compact notation) 



-p(0 )pert 
1 3V 


DDB n 


^P2— 


MlM n B n 


A 2 B r 


(2.17) 


But here the brackets on the l.h.s. are already fixed from a lower stage of the DS problem: the 
self-reproduction conditions for the poles of the gluon self-energy A 2 v(P 2 )- Those enforce, under the 
simplifications adopted, 



p(0 )pert 
L 3V 


DDB n 


} p%=-U r , 2„A 2 


Ur, 2n+1 A 2 1 (jp2 ) j 


(2.18) 


t being the transverse projector of eq. (2A), and u r ^n+i the dimensionless residue parameter in T 
exemplified, for r = 1, by the u\ i3 of eq. (2.5). We conclude that, first, 


[r,0] 

2V 


B [ n ] (P3,Pi) = t{P3+Pi)B\^{p 3 ,Pi) (n = 1, • • .r), (2.19) 


an information more detailed than that of the all-transverse projection fl2.9|) : nonperturbative denom¬ 
inator structure develops only in the variables of transverse gluon legs. Second, 

MlM n = -— tfo). (2.20) 

^r,2n+l 

What has been learned is that as a consequence of the lower (two-point and three-point) equations 
alone , and T^ must contain a term of the form 


- (cfw 


R [r] (^r,2n+l) ^(P) [ r ] 

^ " P 2 +U r ,2nA 2 " 

n =1 ’ 


( 2 . 21 ) 


in the total momentum P, with P 2 = se, of its horizontal channel. But T^y is fully Bose symmetric, 
and therefore must contain analogous terms also for the two crossed channels, with Mandelstam 
variables ue and tp, which we denote by C 2 and C 3 . Thus, 


p[o0] 
1 4V 


-(P'kviv 


(C'Av,2V 


(4'W 


I y[o0] 
+ < / 4V , 


( 2 . 22 ) 


and the derivation shows that the Vfy J defined by this relation contains no more nonperturbative 
denominator structure in se, ue, or tp. The latter conclusion, strictly speaking, follows only for 
the adjoint color representation in each two-body channel, since in Fig. 1 the leftmost gluon line 
projects T^y onto this color subspace. It requires additional considerations, based on the Bethe- 
Salpeter normalization conditions, to check for possible zeroth-order, Mandelstam-variable poles in 
the other colored channels. Since these considerations logically belong to the discussion of the four- 
gluon amplitude, we defer them to the companion article 0. Here we anticipate the result: In 
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zeroth perturbative order , there are no Mandelstam-variable poles in the other color sectors. Thus 

• It Ol 

the ’’reduced” vertex function V^y has nonperturbative structure, rationally approximated, only in 
the variables pj,.. . pf of individual external legs. We emphasize that this in no way precludes the 
existence of glueball-type bound states in color-singlet channels. Such bound states are not elements 
of the generalized Feynman rule T^y at any level r, but arise through the standard mechanism of 
partial ( e.g., ladder ) resummation of quasi-perturbative corrections (g/ 47 r) 2 p rj[^ to all orders p > 1 . 


The poles of (2.15) cannot represent bound states, since in general their residues are not positive 
definite ( at r = 1, for example, ui 3 will be found to be positive ). Thus for the moment they would 
seem to be unphysical artefacts. But one immediately realizes that in fact they play a legitimate 
role by cancelling another unphysical phenomenon. Inspect the analytic structure of the one-gluon- 
reducible terms A t at level [r, 0] (i.e., with [r, 0] diagram elements ). Again suppressing all unnecessary 
arguments, one has 


4 r ’° ] = r^ 0] D [r ' 0] (P) rj^y 1 . 


(2.23) 


By the construction prescriptions for the extended Feynman rules, the gluon propagator contains, 
besides its r + 1 poles, a product of r numerator zeroes of the form 


jPJ (.P 2 + U r ^2n^ 


(2.24) 


n —1 


\ r q1 

On the other hand, both F! S y vertices contain the same product in their denominators, so that there 
remain on the internal gluon line, in addition to the legitimate poles from D describing propagation 
of the exchanged object, a number r of extra poles at positions P 2 = —u rt 2 n A 2 . Such extra poles are 
unacceptable physically; they would imply that the generalized Feynman rule for D is incomplete. 
Now isolate the n-th unphysical-pole piece of A\. From (2.14) and ( 2.1 8| ) , it must involve B n t{P)B n \ 
by computing the residue, the piece is found to be 


B 


fr] (^r,2n+l ) t(P) pH 

n ti 9 , _ * 9 n 1 


P 2 + U rt 2 n A 2 


so that the sum of the unphysical pieces is the negative of ( 2 . 21 ): 

r 

y~^ ( n-th unphysical pole of ) = (Ci^) 2 y 2 v r ' 


(2.25) 


(2.26) 


n =1 


Upon combining (22) and (222) into 


nHo] 

l 4V 


= A' f 


r ’ 0] + Af 0] + A 


/H°] , T/[ r ’°] 

3 “T v 4V ’ 


(2.27) 


the unphysical artefacts then cancel exactly to leave the ’’softened” exchange graphs 

Af° ] = 4’ 0] - (Cl r] ) 2v2v (i = 1, 2, 3). (2.28) 

U Ql 

The poles inferred through the 2U and 3U equations therefore turn out to be ’’compensating 
poles”, cancelling unphysical parts in the one-gluon-reducible terms of T^y. It is clear that in the 
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context of the extended iterative scheme it is the artefact-free Al, rather than the original Ai, that 


represent the physical one-gluon exchange mechanism, and the artefact-free V^y°\ rather than T 


i[r,0] 
4V ’ 


that constitutes the physical extended Feynman rule for four-gluon interaction. On the other hand, 
in the quantity ( 2 . 1 ), 


.,M] = A /[r,0] + A ,[r,0) + v [rfl] _ 


2V,2V ’ 


(2.29) 


an uncompensated C\ remains, enabling, as we have seen, the ( 21)3 term of Fig. 1 to produce the B n 
terms of (|2.14). 


2.3 Extended irreducibility. Rearranged vertex equation 


In the general case where terms (B) 3 , {E) 3 , and (for l > 2 loops) (D) 3 contribute, one first generalizes 
O to the conclusion that the T matrices Tqqw. T' F p vv . T^ v in general have poles (forming 
discrete approximations of cuts) at the same positions, p\ = —u r ^n^ 2 -, in their horizontal channels. 
One then invokes another general property of correlation functions [||: when a pole is present in several 
functions simultaneously, each residue factor is uniquely associated with its own subset of external lines, 
and independent of the remaining legs in the various functions. Thus at the pole p 2 = —u r! 2 nA 2 , the 


residue factor for the rightmost two-gluon configuration in Fig. 1 is the same as in (2.15) for all 
T' amplitudes. Then ( 2.16 ) still follows, but ( |2.17 ) now has four different contributions on its l.h.s., 
corresponding to various dressing mechanisms for the gluon self-energy. Self-reproduction of the B n 
terms with n > 1 in (2.14) is possible if M n and the various matrices taking the place of the Mj in 
(2.17) are all the same mutiple of t(p 2 )- Then the full two-gluon self-consistency conditions [j3j can be 
invoked, which generalize ( 2.1§| ), and the B n term reproduces itself if (2.2C) is imposed. 

As a by-product, one finds that the 1PI functions ^cGWi b ffvv mus t contain terms of the form 


-(dWv 


n[r] (»r,2n+l) ^(-P) R [ r ] 

^ " P 2 +U T)2n A 2 


-( c t') F p,vv = 



(2.30) 

(2.31) 


in their two-body channels (G + G V + V) and (F + F <->• V + V) respectively, where B n , Bp tn 
are the amplitudes analogous to the B n of ( |2.14| ) in the partial-fraction decompositions of the ghost 
vertex and fermion vertices T R (°J R (— p', k,p) with respect to their gluon-leg variable, 

k 2 . The analogous but richer structure in T~ v will be discussed in detail in J|]. 

The structure revealed by these residue-taking operations may look involved at first, but the final 
result is simple: the full 4-gluon, off-shell T matrix (|2.2| ), for example, has no unphysical artefacts 
at all. The artefacts arose because, in a nonperturbative context, the usual decomposition of T by 
the criterion of ordinary one-particle (here, one-gluon) reducibility turns out to be an awkward one: 
both parts in such a division contain unphysical-pole terms that cancel in the sum. It is clearly more 
natural, and better suited to the physics of the problem, to perform the decomposition as in ( |2.27 ), 
where all parts are free of artefacts. To characterize such a decomposition more formally, we call the set 
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of r pole factors common to expressions ( 2.21 ) and ( 2.3C|/fh31 ) a gluonic shadow, described graphically 
by the double wiggly line of Fig. 3, and define as one-shadow-irreducible any amplitude built from 
[r, 0] extended Feynman rules that does not fall into two disconnected pieces upon cutting such a 
shadow line. The defining property of decomposition ( 2.27 ) then is that all its terms are one-shadow 
irreducible. In particular, Vfy J exhibits what one may call extended irreducibility, being irreducible 
both for gluon-propagator poles and for shadow poles, while the ”softened” exchange diagrams A'^’ Q \ 
described graphically by using dotted diagram elements as in Fig. 3, are still reducible for the gluon- 
propagator-poles. 


Use of (2.29) and of the analogous decompositions 

l'h.o] 


j-zi'wj _ y 4 9 L' r > u J _j_ 1 _|_ yMl 


-,/[r,0] _ 

GGVV 
-i/[r,0] 


rpn'^l _ A Ml I | Trh,0] 

1 FFVV ^F ,2 ' ^^,3 ' v FFVV 


GGVV - m GGyV 

- (cfV 


1 FF,VV’ 


(2.32) 

(2.33) 


where A! i and A' F • have obvious meanings as softened exchange diagrams reducible for ghost or fermion 
propagator poles but with ghost-shadow (in non-Landau gauges only) or quark-shadow poles compen¬ 
sated, now leads to the rearranged vertex equation of Fig. 4. With the pole terms of (2.14) having 
been reproduced on the r.h.s. through condition ( 2 . 20 | ) on the four -gluon function, and with the second 

line of (|2.14j) anticipated to be absent by restriction (|2.34|) below, only an equation for the amplitude 
hi 

Bq - an object with three-gluon tensor structure but only two scalar variables - remains. This rather 
strong reduction of the original equation, due to the ” all-in-one-blow” self-reproduction of the B n parts 
with n > 1 through the presence of the compensating poles, generally causes a loss of self-consistency 
conditions and therefore underdetermination, which tends to counteract the overdetermination coming 
from the lack of manifest symmetry. (At r = 1 , for example, the coefficient of ( | 2 . 10 ) appears only 
in and not in B^\ and therefore gets no self-consistency condition of its own). 

In Fig. 4, only the terms that can contribute to the self-reproduction of the extended Feynman rule 
(i.e., produce terms of zeroth perturbative order) at one loop have been made explicit. Thus diagrams 
containing V G q VV or V ffvv no more appear: for these amplitudes, which are superficially convergent, 
it is now true that after extraction of the shadow-pole terms they consist only of superficially conver¬ 
gent loops. We see here that the argument of || concerning these higher amplitudes needs a subtle 
qualification: that argument did not take into account the possibility of certain treelike structures, of 
zeroth perturbative order, which nevertheless appear in the 1PI functions. The shadow-pole terms are 
precisely such structures. Yet their presence does not imply a proliferation of Feynman rules, since 
they consist entirely of building blocks determined already at the level of the basic vertices. For the 
V amplitudes, obtained by subtracting these, the argument of || goes through : their insertion into 
terms (B )3 and (E )3 of Fig. 1 produces integrals for which the number of ^(A^/i/q ) -6 factors lags 
behind the number of Qq prefactors by at least one, and which therefore give only quasi-perturbative 
corrections of order p > 1. Note that this would not have been true for the corresponding T ampli¬ 
tudes before extraction of the nonperturbative shadow pieces. Without this extraction one would have 
missed the dot modifications, and therefore left unphysical artefacts, in the triangle diagrams (A) 3 /, 
(E,E ') 3 ,, and (F,F') 3 , of Fig. 4. 

For the self-consistency problem, the rearrangement for shadow irreducibility brings both simplifi- 
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cations and complications. On the one hand we must now require that the A) of ( [2.28 ), when used as 
diagram-building blocks, preserve perturbative power counting, as did the original Ai by construction. 
It is easy to check from (2.21) that this allows no terms with net positive powers of or p\ in the B n , 
or equivalently, none of the terms in the second line of the p.f. decomposition ( j2.14|) . In the numerator 
polynomial for the complete r.^> rational approximant, we therefore have powers 

C p\) mi {plT 2 {pIP 3 restricted by mi, 2,3 < r, (2.34) 

a restriction significantly stronger than the m ? ; + nrij < 2 r (i 7 ^ j) inferred previously p[], and which 
corresponds to restricting the extra polynomial pj dependence of the spectral functions in (2.12) to at 


most a bilinear one. At r = 1, in particular, this restriction forces the primed coefficients of (2.10) to 


zero: 


x i,2 — 0 


x\ 4 = 0. 


(2.35) 


Moreover it simplifies the writing of the softened one-gluon exchange mechanisms ( 2.28|) , since these 
can now be obtained simply by enumerating the physical gluon-propagator poles with their residues: 
in Landau gauge fixing, 


A 


/[r,0] 


(fci,... k 4 ) = 


J2 [ r 3v(Ah ,k 2 ,p)} 


P 2 = -> 


m =0 


r +iA 2 p? f'" -'—72 [ r 3y(^P, fa, k 4 )] P 2-_ 

Tr,2m+lA pi + a- r)2m+ lA 2 L - 


where 


Pm = [( P 2 + C r r,2m+lA 2 )L> T ( J P)]p2 = _ 


°’r, 2 m +1 


A 2 


<7r,2m+l A 2 5 


rT (^"r,2n ^r,2m+l) 

n= 1 

r ■ 

II (^r,2n+l ^r,2m+l) 

n =0 
(n^m) 


(2.36) 


(2.37) 


and where P = k\ + = — (&£ + k 4 ). In non-Landau gauges, of course, the longitudinal-gluon 

propagator term of the original A\ must be added unchanged. 

On the other hand the rearrangement leads to the result (which is true generally but was not yet 

visible in the simpler case of the two-point equation at one loop treated in Q) that the perturbative 

limit as A —»• 0 cannot, at low levels r, be maintained exactly in all amplitudes, but only asymptotically 

for increasing r. For A —> 0, the v\y^ of ( 2.22| ), as well as the A| r ’ 0 ^ of ( j2.28|) , go over into their zeroth- 

M 

order perturbative counterparts, but the C- for low r do not go to zero: 


[(d rl ) 


1 = _L( 

2V,2VJA=0 p2 \ 


^BM(A = 0)4itL S M(A = 0)} 

Ur-2r,- 1-1 > 


n =1 


Ur,2n+1 


(2.38) 


Thus diagrams involving the one-gluon-exchange mechanisms A), A), A' Fi , such as the triangle dia¬ 
grams of Fig. 4, are expected not to exhibit a fully correct perturbative limit as A —> 0 at low levels r. 
The condition (or rather conditions, because there are several tensor structures involved) for the curly 
bracket in ( 2.3S| ) to vanish, of which we will present examples, can at best be fulfilled asymptotically 
for large r, where they spread their restrictive effect over an increasing number of nonperturbative 
vertex coefficients, and thus become progressively easier to maintain. 
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2.4 3-gluon, r = 1 self-consistency conditions 


To extract the self-reproduction conditions of the extended Feynman rule r!^ one evaluates, with 


eqs. (1.7/1.8) in mind, the divergent parts of the terms on the r.h.s. of Fig. 4 in dimensional regulariza¬ 
tion and Landau gauge fixing, with [1,0] diagram elements throughout. The restriction to divergent 
parts makes these calculations somewhat analogous to the computation of one-loop renormalization 
constants in perturbation theory (and much more feasible than full evaluations of [ 1 , 1 ) radiative cor¬ 
rections, which already at r = 1 are very lengthy). For the input to diagram (C) 3 /, we anticipate 
formulas from appendix B of the companion article j|]: this approximant represents a theoretically 
motivated restriction to a subset of fifteen of the many possible color and Lorentz tensor structures of 
a four-gluon amplitude, with invariant functions characterized by a set £ = {( 1 ,... £ 17 } of seventeen 
dimensionless, real numerator coefficients. The terms proportional to the number Np of quark flavors, 
arising from fermion-loop diagrams (F) 3 / and (F') 3 i, are valid for massless quarks (mjr = 0 in the 
notation of the appendix of Q), where all fermionic mass scales, too, are simply multiples of A. For 
brevity, we abstain from listing contributions of the various diagrams separately [jEj] and present only 
the combined results. Eqs. ( 1.1C| ) for the coefficients xqj (now written Xi for brevity) of the Fq 1,0 ^ 
invariant functions (| 2 . 10 |) are: 


1 

Ao - 


9 15 1 (\ \ 

- - 7 %i + —X 3 H- -xix 2 - 9xix 4 + x 3 x 4 

4 16 u 3 V4 / 

2 / 1 

+ P N p[Z 3 “I- X 4 Z 3 

3 \ u-x 


-2424 

w 3 


= X 1 


(2.39) 


1 I" 3 

A) 1-2 


7*3 


1 (L n It> t> \ 

H-77^2^4 - 2x 4 - —X1X5 + 7X3X5 

u 3 \2 2 4 / 


15 


+ o N f[4 H-x 5 z 3 - 24 

3 V u 3 w 3 


= X 2 


(2.40) 


1 

A) 


r 3 

1 ( 

37 

3 \ 

77X3 + 

— 

--X1X4 + -X3X4 

12 

u 3 \ 

4 

2 / 


2 A7 / 1 1 \ 

+ - N F z 3 H-x 4 2 3 - 2424 

3 V u 3 w 3 / 


= x 3 


(2.41) 


1 

Po 


r 9 

15 

1 / 

31 

5 5 \ 

- 7*1 + 

77*3 + 

— 

--X1X2 - 

7x1x4 + 7x2x3 

L 4 

16 

u 3 V 

4 

4 4 / 


2 / 1 1 

+ —Np 23 H- x 2 2 3 - 2424 

3 \ u 3 w 3 


= X\ 


(2.42) 


1 


1 


1 


15 


A L2' ,: l + V 3 ( - i X2Xi - T* “ + l X3Xs ) - Zl(<) 

+ + —X5Z3 - —zl 

3 V u 3 w 3 


= x 4 


(2.43) 


(The fermionic vertex-coefficients 2 * are defined by eq. (L 6 ) below. 
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On the other hand, the terms of the invariant function F\ of ( 2.11 ) turn out to be fed only by 
themselves, and by terms in Viv which our above-mentioned, restricted form of that amplitude omits 
in the first place. Again, since we view the basic vertices as an interrelated whole, it did not seem 
consistent to us to keep just one source of such terms. Thus 


x 6 = 0 


xy = 0 


(2.44) 


is a consistent and self-consistent choice in our framework. 

A noteworthy feature is that the coupling to the 4-gluon amplitude enters only into the three 
equations ( |2.41[ - |2.43|) , and only through two linear combinations of its seventeen coefficients C, 


Si(C) = ^(3Ci-Cr), 

Z 2{() = — ( 3(2 + 3C3 — C8 - C9 ) • 


(2.45) 

(2.46) 


The observation that the rather large (and, as it will turn out, strongly overdetermined) self-consistency 
problem of the four-point vertex parameters couples to the 2-point and 3-point problem only through 
this narrow ’’bottleneck” will be important as it will suggest ways of breaking down the rather volu¬ 
minous total self-consistency problem into more manageable pieces. 

As already noted, there is no equation with X 5 on its r.h.s. But still has three-gluon tensor 
structure, and xi, by (fO|), appears twice in conjunction with two different tensor structures, so there 
are two equations, ( 2.39| ) and ( 2.42 ), for x\. The relation obtained by subtracting these, 

— ( 8 x 1 x 2 - ^-xix 4 + X 3 X 4 - 7 ^ 23 : 3 ) + Zi(() - 7 N f (— — —\z 3 = 0, (2.47) 

u 3 \ 4 4/ 3 V u 3 / 

represents the imposition, in zeroth perturbative order, of Bose symmetry on a DS equation that is 

not manifestly Bose symmetric. It appears to be fortuitous that the ’’loss” of one equation, incurred 

in the reduction of the self-consistency problem to the partial amplitude is just compensated by 

one Bose-symmetry restriction; we are not aware of a deeper reason for this phenomenon. 

Finally we note the result for the coefficient of the perturbative-remainder divergence, the of 
eq. flOl ): 


~^3 T ) — 0 


(,P2 ~P3) P 

+6 ap ( P3 - P1 r 

+ 5 pK ( Pl - P2 r 


17 1 

4 u 3 
17 J_ 
4 u 3 
17 1 

4 u 3 


~ x, - 8 x 1 x 3 + 7 X 3 + -N f 1 H- x 3 z 3 - z 1 

4 / 3 V u 3 w 3 


37 2 3 

- ~ X 1 + ^x r x 3 

37 


1 2 

— 4 

w 3 

- -jr x \ + 7X1X3 ) + 7 N f ( 1 + — x\z 3 - —z\ 
4 2 / 3 V u 3 w 3 


2/1 

+ -1Vf(1 H- x±z 3 

3 V u 3 

1 


The purely perturbative result would be of the form (1.11) with the Landau-gauge value 


ii), 


17 


‘3y(^ - 0) - -~£ + % Nf - 


(2.48) 


(2.49) 


The existence, and lack of complete Bose symmetry, of deviations from the perfect perturbative limit 
do not come as a surprise in view of what we noted before in connection with ( [2.38| ). Eliminating 
them would again produce overdetermination of the x coefficients and generally is not feasible exactly 
for low r but only asymptotically for large r. 
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3 Equations for the ghost vertices 


For the ghost-gluon-antighost vertex, Tq V q, and its generalized Feynman rule, it is again self-consistent 
(though not the most general solution) to assume a pure / a & c color structure, the Lorentz structure 
then being given by 

[rj$j(p, k, -p')Y abc = if abc [p^F^ 0] (p, k, -p') + k^F[ r ’ 0] (p, k, -p')\ (3.1) 

The dimensionless invariant function F), with perturbative limits p'j°' ,pert _ depend on the in¬ 
variants p 2 , k 2 , p /2 . Fig. 5 shows the diagrammatic form of the DS equation for T gvg in its ghost 
channel, i.e., with the ”G” leg as the unsymmetrically distinguished leftmost leg. It again features a 
four-point amplitude T' (:vv( - ; which is 1PI in only the ’’horizontal” channel. Residue-taking both in 
this equation and in the corresponding equation in the antighost channel again reveals the presence of 
compensating poles in T\ and taking these into account one again obtains a rearranged form of the 
equation as in Fig. 5(b). 

When staying strictly in Landau gauge, as we do in this paper, it is actually unnecessary, as far as 
the self-consistency of the generalized Feynman rule is concerned, to evaluate the terms on the r.h.s. 
of Fig. 5(b) in detail: brief inspection shows that the latter, at £ = 0, do not sustain nonperturbative 
A terms. Consider e.g., term ( B)q of Fig. 5(b) with the momentum assignments shown. Its upper 
^ GVG vertex has, by (|3d|), Lorentz structure 

ip' - q- 2 ) X F 0 + q$F 1 . (3.2) 


The internal gluon line carrying momentum g 2 has, in Landau gauge, a transverse projector f KA ((? 2 ) 
so that only the p ,A Fo portion survives. But p’ is an external momentum not running in the loop, so 
the integrand loses one power of the loop momentum as compared to standard power counting. Since 
the loop had only logarithmic divergence to begin with, it is now actually convergent. Analogous 
arguments apply to the term ( C)q of Fig. 5(b). The only possible nonperturbative modifications are 
therefore those from term ( A)q , if any. These must have at least one denominator factor (p 2 + 12) 2 A 2 ) 
in the momentum variable of the leftmost (ghost) leg. But the same argument applies to the two 
alternative forms of the ghost-vertex equation not displayed in Fig. 5, with the ’’antighost” (momentum 
— p') and gluon (momentum k ) lines, respectively, as leftmost legs. The only nonperturbative terms 
in the F invariant functions of equation ( |3.l|) , that are candidates for self-consistency, are therefore 
those proportional to 

_ A? _ 

(p 2 + u[ 2 A 2 ) ( k 2 + u[ 2 A 2 ) (p /2 + u[ 2 A 2 ) 


These terms, however, make the loop appearing in term ( A)g of Fig. 5(b) convergent. Since the self¬ 
reproduction mechanism of eqs. (|1 .7/1.8) is dependent upon the j divergence factor, the amplitude 
cannot develop any zeroth-order nonperturbative terms in Landau gauge: 


[ ] r GVQ(P , ~P')] abc = ^fobcpP = (r aveTabc ^ = °> a11 ^' 


(3.3) 
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The effect is, of course, basically familiar from perturbation theory: there, the one-loop divergence 
of the renormalization constant Z\ vanishes at £ = 0. The preceding discussion merely serves as a 
reminder that in the present context such special divergence reductions also have qualitative dynamical 
consequences as they suppresss the divergence-related self-consistency mechanism. For an amplitude 
with unphysical degrees of freedom such as T GV q, it is of course legitimate to depend on the gauge 
fixing in this way. 

The ghost-self-energy equation, due to its general divergence reduction as discussed in appendix 
A.2 of also has effectively a logarithmically divergent integral. When evaluated with the purely 
perturbative vertex ( |3.3| ), that integral can produce no more than the perturbative divergence, so that 
again no nonperturbative terms are formed: 

D^(p 2 ) = D^ pert (p 2 ) = \ (£ = 0, all r). (3.4) 

pZ 

We see that the assumption of refs. [j5j that ghost vertices remain perturbative is justified only in 
Landau gauge. 


4 Massless-fermion vertices 


In the absence of Lagrangian mass terms for quarks, nonperturbative mass scales in the two basic 
fermion vertices T FF and T FVF can only be multiples of the A scale. This case is technically far 
simpler than the rather complicated situation encountered in the presence of additional RG-invariant 
scales from ’’current” quark masses, and is the only one we consider in this paper. 

The DS equation for the inverse quark propagator is given diagrammatically in Fig. 6. The 
corresponding extended Feynman rule at level r = 1, as discussed in the appendix of Q, is 

<4 ' 1} 

Its nonperturbative content is characterized by the three dimensionless, real parameters w \, w 2 , W 3 . 
Extraction of the self-consistency conditions for these involves calculating the divergent parts of the 
loop of Fig. 6, evaluated with [1,0] input elements, and proceeds largely as in the gluon-propagator 
case considered in Q. The resulting equations 


1 

A> 

1 

A) - 


4wi — Az\ 
4it>i2:i — 4^2 


w 2 = w 2 , 

= , 


= W3 , 


(4.2) 

(4.3) 

(4.4) 


make reference to the parameters 2 of the quark-gluon three-point vertex, T FVF . Its transverse-gluon 
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projection at level r = 1 reads, 


[r ££L(-p'>Mr= t^( k ) + z i + A W / a 7 m + 7 a 


A 


ji + u^A 


+ "V + ^A 7 ^ + w'A + k 2 + u' 2 A 2 + Z 4 (jl' + ^A 7<< + 7 % + A 


+ ^5- 


A 


A 


-sy r 

’ f + w 2 A f> + w 2 A. 

Here the notation of ref. ||] for the dimensionless coefficients has been changed and simplified some¬ 
what, the relation | 3 ] —> this paper being given by 


(4.5) 


jn 

z 0,l 


zi, 


7 w 

0,4 


* 2 , 


,[!] 

- 1,0 


*3, 


7 W 

- 1.1 


Z4, 


7 W 

h,4 


Z5- 


(4.6) 


We have from the outset omitted all terms that would lead to conflict with perturbative divergence 
degrees. For T FVF there are two DS equations, one in the ”fermionic” and one in the ’’gluonic” channel, 
which are equivalent for the exact vertex but in general will give rise to different approximations. 
Rearranged for one-quark-shadow irreducibility in the now familiar way, these are depicted in Figs. 
7(a) and 7(b), respectively. While each of the two forms, due to the compensating-poles mechanism, 
suffers from ’’loss of equations” in the sense discussed in sect. 2.3, it is interesting that the two forms 
taken together produce just the required number of self-consistency conditions: 
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In addition, the quasi-perturbative remainders contain divergences given by 
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= z 3 


= Z 4 
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for Figs. 7(a) and 7(b) respectively, which differ from the perturbative quantity 

4vf(£ = 0) = \ 


(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 


by defect terms involving the vertex constants x and z, which again cannot be forced to zero at r = 1 
but only asymptotically at large r. 

Note that there is no equation determining the vertex-pole position w' 2 , which by (4.2) is also the 
propagator-zero position. 
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5 Solution for 2- and 3-point coefficients 

5.1 Analysis of equations 

The system of self-consistency conditions for the r = 1 nonperturbative coefficients, as established up 
to now, consists of eqs. (4.8/4.30/4.31) of ref. [!(| for the gluonic self-energy parameters ui,u 2 ,u 3 , 
which now assume^ the simpler forms 

u' 2 = u 2 = u 2 , (5.1) 


1 

A) 


9 


33 


1 


-Ul -— X\ + -X 3 - 2 Nf W 3 + (wi + W 2 )z\ + Z 2 - ~ Z S 


U\ , 


(5.2) 


1 

Po 


-U 2 


(3xi - ^x 3 ) + 9uix 3 - 9x 4 - 2N f (z 3 (w3 + ^u 2 ) + (w 2 - Wi)z 4 + 2:5 


= u 3 


(5.3) 


plus eqs. ( 2.39 - T43| ) above for the 3-gluon-coefficients X 1 ...X 5 , plus eqs. (4^-4. 4) and (4.7- 4.11 ) above 
for the self-energy coefficients W 1 ...W 3 and vertex coefficients z\. ..25 of massless fermions. Its peculiar 
properties, in particular with respect to under- and overdetermination tendencies, can be summarized 
as follows. 

(a) One observes that the system as a whole exhibits a scaling property: any one of the coefficients 
that is presumed to be nonzero may be divided out of these equations, while replacing the others by 
their ratios to this one or to a uniquely fixed power of it, and rescaling A accordingly. This property, 
which in a nonlinear system is nontrivial, is a natural consequence of the scheme-blindness of the 
basic self-consistency mechanism: a rescaling of A, which corresponds to a change of scheme, will not 
change the form of the zeroth-order conditions. We will choose, for definiteness, a rescaling by the 
3-gluon coefficient xi of fl2.10|): 


A 2 = xiA 2 , 
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Zl = 

Zl 
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z 2 

Z3 = 

Z3 

= 

ZA 

Z5 = 
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1 

Xl 

X\ 


x\ 


(5.4) 


(5.5) 


The reason for this choice is that by putting xi = 0 one would end up with only the trivial solution (all 
nonperturbative coefficients vanishing), so one is not losing interesting solutions by assuming xi A 0. 

x Eq. ( |5.2| ) corrects for a misprint in eq. (4.31) of the first of refs. Q, where a tadpole contribution —97Vcui,i/4 
appears with the wrong sign 
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( For the scaling of the fermionic parameters, we are for the moment assuming x\ to be positive ). In 
effect, the rescaling reduces the number of unknowns by one while introducing the modified A of (|5.4|), 
scaled by an unknown factor. 

The rescaled system of fourteen conditions, with fifteen unknowns and Z i, Z 2 as external param¬ 
eters, has solutions coming in pairs: one checks that if 


{ 


Ul,U 2 ,U 3 , W ll W 2 ,W 3 , X 2 ...X 5 , Z 1 ,Z 2 ,Z 3 ,Z 4 ,Z 5 


is a solution, then for the same Z\ and Z 2 the set 

{ 


Ul,U 2 ,U 3 , -W^-W^Ws,, x 2 ...x 5 , -z 1 ,z 2 ,z 3 ,-z 4 ,z 5 


(5.6) 


(5.7) 


is also a solution, which we shall refer to as a ’’mirror” solution. Note that this discrete ambiguity 
affects only fermionic parameters. 

(b) We have not obtained equations fixing the vertex-denominator parameters (u 2 and w 2 in the 
present case). Neither analysis of the 4-gluon vertex ||] nor, as preliminary studies indicate, use of 
’’resummed” DS equations will change this situation. We did obtain conditions like the w 2 = w 2 of 
eq. (4.2), and the corresponding u' 2 = u 2 = u 2 of eq. o, which ensure one common pole position 
in all basic vertices for a given type of external leg (and also the presence of propagator zeroes at 
the positions of vertex poles), but u 2 and w 2 . in the end, have no determining equations of their 
own. This leads to the unexpected conclusion that the divergent parts of the momentum-space DS 
equations as used up to now do not yet determine the nonperturbative A dependence completely. The 
reason is that these equations, in a sense, do not provide enough divergence. Indeed, the quadratically 
divergent gluon self-energy is the only vertex having u 2 and w 2 appear at least on the right-hand sides 
of its self-consistency conditions for u\ and u 3 , but at least two more equations with the same or 
higher degree of divergence would be needed to “lift” the two parameters from the denominators of 
loop integrands into numerator expressions that provide self-consistency conditions - a feat that only 
divergent integrations can perform. 

Additional conditions for fixing u 2 and w 2 therefore should have general compatibility with the 
momentum-space DS equations and provide sufficient divergence. The only natural candidates here 
are those requiring the vanishing of the ” equation-of-motion condensates”, i.e. of vacuum expectations 
of the simplest (dimension four) local composite operators proportional to the left-hand sides of the 
field equations. We use the condensate conditions for the ghost and fermion fields in the form 


(90^0) <0| C a (x) {[dafeD + goVofabcAf (x)d»] C b (x)} |0) = 0 

( fi ' ot ' o ) 2 (°l $( x ) {[*0 + 9oVo 4( x )\ H x )} 1 °) = 0 


(5.8) 

(5.9) 


In momentum space these are, of course, nothing but the ghost and quark propagator equations in¬ 
tegrated over momentum space, or equivalently, taken at zero separation in coordinate space. They 
are therefore obviously compatible with, and natural completions of, the unintegrated ( momentum 
space ) or uonzero-separation ( coordinate space ) DS equations we have exploited up to now. In stan¬ 
dard integral-equation theory with convergent integrals and well-behaved functions, they would not 
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represent independent statements, but in a theory with divergent loop integrals and therefore in need 
of renormalization, they do carry new information: since they involve operators of higher compos¬ 
iteness, they possess new divergences leading to new zeroth-order conditions on the nonperturbative 
coefficients. At the r = 1 level these conditions read, 

4 [-«? + u 3 ] A 4 = 0 (5.10) 

Po 


_3 

Po - 
2 


wf — Awfw 3 — 2 w\W 2 W 3 — w\w 3 + VJ$ 


O A O 

U 3 — u l — 3uiw 3 — 6 w 1 + 12w 1 w 3 + 6 w\W 2 W 3 


+ 2>zi{u\W\ — U 1 W 2 + 3w\ + 2w\w2 + 2w\w\ — Qw 3 w 3 — 2w2W 3 ) 


+ 3z 2 (—iti - 2u;f - 2w\W 2 — 2w\ + 2w 3 ) + z 3 (u\ + u 2 + 3w 3 ) + 3zpW2 - w\) + 3 z 5 


A 4 = 0 


(5.11) 


Condition ( |5.10| ), from the ghost equation of motion, notably provides a restriction on gluonic 
parameters ( it is, incidentally, equivalent to requiring the vanishing of the zeroth-order, dimension- 
two gluon condensate (A^A^), which at r = 1 turns out to be proportional to —u± + u 3 ). Since 
in Landau gauge at one loop it is the only condition from the ghost sector, is independent of the 
presence of fermions, and of remarkable simplicity, we give it priority in complementing eqs.(A2|/5.3). 
Condition (5.11) brings in, in addition, the fermionic parameters, and is suitable for complementing 
eqs. (4.3/4.4). The order-[l,0] equation-of-motion condensate for the gluon field, 

(So^) 2 <0| [l(a<x - a" as ) 2 + 2 ( 9 „„5) - a"At) A£At 

+ (go^of fabefcdeKK A[l c A v d + (g 0 l/%) /abc(<9 M C 0 ) A%C b + (fi^o) N F p4P\ |0) = 0, (5.12) 

is the most complicated, and its zeroth-order form at r = 1 , 


r 24 
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u x - U 3 
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N f 
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+ 


— 390uj + 228u 3 — 324x4 + 594uiXi + 279uix 3 + 270u 2 xi — 45 u 2 x 3 


8uf — 8u 3 + 24uiu; 3 + 48 w\ — 96 w\w 3 — A&W\W 2 W 3 


+ 24 ^i (— u\W\ + u\W2 — 4 u>f — 2w\w2 — 2w\W2 + 6 w 3 w 3 + 2w2W 3 ) 

+ 24z 2 (ui + 2w\ + 2w\W 2 + 2w% — 2w 3 ) — 8z 3 (ui + u 2 + 3 w 3 ) + 24 za(wi — rc 2 ) — 24^5 

1 

+ M 

involves the largest number of vertex parameters simultaneously. In the present context it cannot be 
applied in its exact form; the 3-loop-terms must be omitted for formal consistency with the omission 


3-loop-terms 


A 4 = 0, (5.13) 


21 






















of 2-loop-terms in the l = 1 gluon-self-energy calculation. For these reasons, condition ( 5.13 ) may be 
expected to be the most difficult to fulfill on our level of approximation, and is not one of our primary 
choices for completing the self-consistency system. We will check in the end to what extent it can be 
accommodated. 

(c) Due to the large dimensions and considerable overdetermination of the 4-gluon self-consistency 
problem to be discussed in Q, the total coupled problem cannot, in our experience, be attacked directly 
with currently existent mathematical software tools. However, we have already emphasized the ( also 
unexpected ) result that the 4-gluon-vertex problem couples to the fewer-point amplitudes only through 
the narrow ’’bottleneck” of two 4-gluon-coefficient combinations ( |2 . 45| / pT46|) appearing in only three 
of the 3-gluon conditions. This situation of a near decoupling of the 4-gluon self-consistency problem 
renders the following strategy sensible ( it is, in any case, the only practical strategy at present ). One 
omits, as a first step, the 4-gluon conditions completely, and treats the parameters Zi, Z 2 appearing 
in ( |2. 41 -2.43) as two additional unknowns in the 2-point-plus-3-point-system, which thereby becomes 
doubly underdetermined. Combined with the scaling property of point (a), which effectively reduces 
the number of unknowns by one, this results in an effective one-parameter freedom in the solution of 
the 2-plus-3-point problem. Since the number of calculable coefficients - in the present case, fourteen 
in the 2-and-3-point amplitudes, and seventeen in the ’’minimal” four-gluon vertex to be discussed in 
[j|] - is much larger, the solutions will still be nontrivial and informative; in particular, one may explore 
in what range, if any, of this one-parameter freedom there exist physically acceptable solutions. 

In a second step, which we defer to H, one may then adjoin the values of Z\, Zo thus determined 
as additional constraints to the 4-gluon self-consistency problem: this will represent only a minor 
increase in the anyway massive overdetermination of that problem. Since the 4-gluon system refers 
to the 2-and-3-point coefficients, it inherits the effective one-parameter freedom, and to within that 
freedom may be dealt with separately, with methods adapted to its overdetermined nature. 

It would seem that any one coefficient or combination of coefficients of the 2-and-3-point system 
could be used to parametrize the effective one-parameter freedom; in particular, some combination of 
the quantities ( 2.45 / 2.41 ), which caused the freedom in the first place, would seem to be a natural 
parameter. However, one again faces unexpected restrictions here: due to the peculiar structure of the 
system in its fermionic unknowns, the fixing of a combination of non-propagator parameters, instead 
of rendering the system well-determined, usually splits it into an over- and an underdetermined part. 
The parametrizing quantity should thus refer to propagator coefficients. In the following we will 
choose, for no other reasons than technical simplicity, the rescaled quark self-energy coefficient w\. 

(d) At present, the question remains open as to whether there exists a preferred or natural way of 
finally removing the one-parameter freedom. One might think of recalculating the quantities Z\, Z 2 
later from the least-squares four-gluon solution and see if there is a parameter range where they agree, 
at least qualitatively, with those from the 2-and-3-point solution. We shall indeed do this in J|], but 
shall see that in the parameter range where the entire solution is physically acceptable, a mismatch is 
unavoidable at r = 1, although small in the case of a pure gluon theory. Alternatively, the vanishing 
of any of the previously noted approximation errors existing at the (r = 1, l = 1) level could be used 
as a condition. The common problem of all conditions of this kind is that (i) there are several of 
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them, and any selection from among them appears arbitrary, (ii) they are mostly so restrictive that 
their imposition leaves only the trivial solution, with all nonperturbative coefficients vanishing. The 
message the defect terms seem to convey is that the still rather simple and rigid structure of the 
r = 1 system of approximants entails unavoidable approximation errors that cannot be forced to zero 
without overstraining that structure; they can disappear only gradually as r is increased. 


5.2 Discussion of solutions 


The system augmented by ( [5.10 ) at a fixed value of w\ may be reduced by successive elimination 
to an algebraic equation of the 10th degree for the quantity W 3 . The other coefficients can then be 


calculated recursively from the solutions of this equation and eq. ( 5.11 ), and depend parametrically 
oniii. ( These calculations have been performed using the MAPLE V computer-algebra system ). 
The following noteworthy features emerge. 


(a) When assuming x\ < 0 and performing the rescaling (545/ 5A) with |aq| = —x± instead of 
xi, the ten roots obtained for W 3 are all complex. Such solutions can immediately be discarded as 
unphysical, since they lead to vertex functions not real at real Euclidean momenta, and the nonlinear 
nature of the system permits no superposition to obtain real solutions. Therefore no physical solutions 


have been lost by assuming x\ > 0 and rescaling as in (5.5/ 5.4). 

(b) Over a range 0.3 < w\ < 1.2 ( all ranges quoted are approximate ) only eight of the ten w 3 
roots come in complex-conjugate pairs, but two are real: there exist solutions with all vertex coefficients 
real. This result is entirely nontrivial, and represents substantial progress over the earlier attempt of 
refs. |J, where what we would now call the r = 1 level of approximation was studied in a more 
heuristic fashion, with strong a priori simplifications of the vertex approximants, and without taking 
the compensating-poles mechanism into account. There, only partly real solutions could be found. 

Of the two real W 3 roots, one is negative and one positive. The negative W 3 value always turns out 
to lead to ’’tachyonic” pole positions (negative values of the p± of ( 5.1 8|) below) in at least one of the 
two propagators, and can also be discarded as unphysical. Again it is nontrivial and noteworthy that 
only one of several solutions of the nonlinear system stands out as a candidate for a physical solution. 
Each of the two real solutions still exhibits the doubling of eqs. (5.6/|5.7), i.e. has a mirror solution 
for some of its fermionic parameters in the range —1.2 < w\ < —0.3. 

(c) For the solution in the range 0.3 < w\ < 1.2 with W 3 real and positive, table 5.1 records the 
nature of the poles in the Euclidean transverse-gluon and fermion propagators, which now read 


D [ l'°\k 2 ) = 


k 2 + 142 A 2 


( k 2 + a + A 2 )(k 2 + cj_A 2 ) ’ 

S M«) =_ * + "* A _ 

( see eqs. (| 2 T|) and ( [PI ) ).One finds that over a narrower range of w\, namely, 


(5.14) 

(5.15) 


0.5 <m< 0.9 


(5.16) 
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0.3 < n)i < 0.4 

0.4 < w\ < 0.5 

0.5 <w\< 0.9 

0.9 < w\ < 1.1 

1.1 <wi< 1.2 


< 0 

complex 

< 0 

< 0 

d_ 

> 0 

conjugate 

< 0 

< 0 

fir 

> 0 

> 0 

complex 

> 0 

~P 2 - 

> 0 

> 0 

conjugate 

> 0 


Table 5.1: Ranges of interest for propagator-pole parameters 


the solution with the real and positive w 3 root in addition fulfills the three inqualities 


7v = u 3 


U\ - U '2 


> 0 , 


7 F = w 3 


W 1 — W 2 


> 0 , 


\w\ + w 2 \ > 0, 


(5.17) 


so that both propagators simultaneously exhibit complex-conjugate pole pairs at k 2 = — <r±A 2 and 
p 2 = —p\ A 2 , where 


u\+u 2 

<?± =-X-± *7V7 


p± = 


w\ + W 2 


— w 3 ± i(w 1 + w 2 )'tF- 


(5.18) 


Thus there exist solutions in which the elementary excitations of the two basic QCD fields are both 
short-lived. In the present framework this is the essential indicator of confinement, since it implies the 
vanishing of S-matrix elements with external single-gluon or single-quark legs J5|. We again regard 
it as nontrivial that a parameter range should at all exist in which this situation prevails. Note also 
that in the gluonic portion of table 5.1, there is always at least one “tachyonic” gluon-propagator pole 
outside the slightly wider range 0.4 < u>i < 0.9 ); we view it as significant that the only solutions with 
real vertices and non-tachyonic gluons have gluon propagators with complex pole pairs. - Over the 
range (5.1 6|) , all other vertex coefficients are only weakly varying. 


This interesting solution still has a ’’mirror” solution in the sense of (5.7), i. e. in a w\ range which 
is the negative of ( [5 ■ 1 6| ) . We are not aware of a theoretical criterion that would resolve this discrete 
ambiguity. One might prefer the solution with negative w\ on the empirical grounds that it gives the 
r = 1, zeroth-order fermion condensate 


[b°] 12 , o \ n 

= — (w 1 - 2wiw 3 - 2w 2 w 3 ) A 
Po 


(5.19) 


the negative sign established in the context of current algebra und QCD sum rules. However, there is 
no reason for believing that a solution with low r, which is generally crude and more so in the fermion 
sector, must already give the correct sign for such a sensitive quantity. 


(d) Upon imposing condition (5.13), from the gluonic equation of motion, to remove the residual 


freedom in the w\, we find values w\ not only outside the range (5.11) but in fact outside the larger 
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range of table 5.1 , where w 3 and thus the entire solution creases to be real and physically acceptable: 
enforcing ( |5.13| ) one loses the possibility of a physical solution and of all the features noted in (c). 
Thus ( |5.13| ), like removal of the order-g 2 defect terms noted above, seems to be a strongly restrictive 
condition that the simple r = 1 structure is too rigid to accommodate. 

Within the general strategy suggested and used here, it appears that restriction to the quite 
limited parameter range where all propagator singularites are complex conjugate, and none tachyonic, 
in itself represents a sensible limitation to the one-parameter freedom, and one that is difficult to 
narrow further without overburdening the r = 1 approximation. 


Ui U 2 U3 

Xi X2 X 3 X4 X 5 

iv 

-0.3604 -0.4884 +0.1299 

+1.0000 -8.7433 +8.9088 -3.2607 -6.2711 

0.3547 

W\ w 2 w 3 

Zi Z 2 Z 3 Z4 z 5 

7 F 

+0.6749 +0.6749 +0.1202 

-0.9561 -0.9356 -0.4282 +0.4094 +0.2242 

0.3468 


Table 5.2: Typical solution with Np = 2 


(e) For use of the generalized Feynman rules in applications, we list in table 5.2 a typical set of 
two-point and three-point vertex coefficients for Np = 2, for w\ chosen in about the middle of the 
range (|5.16 ). ( Since in that range W 2 varies slowly and is itself of modulus ~ 0.7, we choose the 
point w\ = W 2 for simplicity.) This set still needs completion through the corresponding four-gluon 
vertex coefficients but since many lower-order calculations need at most 3-point vertices, it seems 
legitimate to defer presentation of these to |1. The propagator-pole parameters (5.IS) for this solution 


are 


a± = (-0.4245 ± ffi.3547)A 2 


pi = (0.3353 ± ffi.4679)A 2 


(5.20) 


For purposes of comparison, we also briefly look at solutions for the pure-gluon theory (Np = 0). 
Here the parametrizing quantity may be taken to be the gluonic vertex coefficient x 3 , and again we 
choose a typical value, £3 « 1 , from the ( again existing ) range in which the gluonic propagator poles 
are complex conjugate. ( The value x 3 = 1, incidentally, is also one which symmetrizes, though not 
removes, the defects of ( [2.48 ) in the perturbative three-gluon divergence ). Table 5.3 lists coefficients 
for this case. It is still impossible here to accommodate condition (|5.13) in a physically acceptable 


U\ U 2 u 3 

Xi X2 x 3 X4 X 5 

lv 

-1.7429 +0.8456 +3.0376 

+1.0000 -6.1825 +1.0000 -4.8682 +28.605 

1.1650 


Table 5.3: Typical solution for pure-gluon system 


solution. Note however that this solution is now 


unaffected by the doubling of & 
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The presence or absence of the massless-quark loops obviously has a strong effect on several of the 
coefficients. In particular, the value of the transverse-gluon propagator function ( |5.14| ) at k 2 = 0 - a 
finite constant whose sign is determined by the 112 parameter - is positive for the pure-gluon system 
but negative in the presence of the light quarks. It is again unlikely that the crude r = 1 solution 
should describe this effect quantitatively, but its qualitative trend is plausible from the minus signs of 
the fermionic self-energy loops. 

6 Conclusion 

We have demonstrated the feasibility of a self- consistent determination of generalized Feynman rules, 
accounting for the nonperturbative A dependence of correlation functions through a modified iterative 
solution, at the simplest level of systematic approximation of that dependence. We have shown, 
and regard it as nontrivial, that the nonlinear self-consistency problem admits physically acceptable 
solutions, that these stand out clearly against a majority of unphysical ones, and that there exist 
solutions in which both of the elementary excitations of the basic QCD fields exhibit the short-ranged 
propagation described by complex-conjugate propagator poles. 

It is useful to recall the restrictions under which we have studied this self-consistency problem. We 
have considered the r = 1 level of rational approximation of the A dependence - the lowest level of 
interest for a ’’confining” theory like QCD. The limitations inherent in this low approximation order 
have become clearly visible; its structure is far too simple and rigid to satisfy all desirable conditions 
and restrictions simultaneously. 

We have worked with the ’’ordinary” DS equations only, with bare vertices on their distinguished, 
left-hand external legs. We have evaluated the self-consistency conditions on the one-loop level, in 
Landau gauge, and with a special decoupling of the 4-gluon-conditions as suggested by the peculiar 
’’bottleneck” structure of the system. It is desirable for future work to gradually remove these limita¬ 
tions, and in particular to study the Bethe-Salpeter-resummed forms of the vertex equations, which 
may have more of the important physical effects shifted into the low loop orders. 

Even with such improvements, two problems are certain to persist that have emerged clearly from 
the present study. One, which arises only when studying vertices with at least three legs, and has 
nothing to do with the specifics of the present method, is the overdetermination dilemma unavoidable 
when seeking approximate-but-symmetric solutions to the not manifestly symmetric DS equations. 
This may be “swept under the rug” by trivial symmetrizations, but only at the expense of depriving 
oneself of an important measure of error. The second problem is that the DS equations, through their 
divergent parts, do not fix the common set of denominator parameters of the approximants. A way 
of understanding this interesting result is to recall the relation with the operator-product expansion, 
as discussed in sect. (2.3) of ||: the OPE, in its higher orders, contains vacuum expectations of local 
operators of arbitrarily high compositeness, whereas the DS equations contain at most insertions of 
three operators at the same spacetime point. The extra composite-operator renormalizations required 
by those higher condensates represent extra information which the usual DS system does not supply 
directly. Equation-of-motion condensate conditions, which do represent statements about quantities 
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of higher compositeness, are capable of supplying the extra information, and natural complements 
insofar as they are special, zero-separation cases of DS equations. Their role, which in the present 
context may still have looked marginal, will clearly become more central when going to higher levels 

r. 

It should be kept in mind that for a confining field system such as QCD, the generalized Feyn¬ 
man rules as considered here allow only the calculation of off-shell Green’s functions of the elemen¬ 
tary fields. These still carry little observable information, although the spectrum of the elementary 
excitations as determined by the singularities of the two-point-functions does constitute important 
qualitative information. To calculate on-shell amplitudes, whose external legs are bound states, one 
would need in addition bound-state vertices to sit at the outer corners of S-matrix diagrams. These 
have not been touched upon in this paper, since they are conceptually quite different from the zeroth- 
perturbative-order quantities: they arise from partial (ladder or improved-ladder) resummation of 
quasi-perturbative corrections g 2p V^’ p \ with p > 1, for certain superficially convergent amplitudes 
Tjv, in which the mechanism of eq. ( |1.8| ) plays no role. Their determination must therefore rely on 
the established Bethe-Salpeter methods for bound states. 

On the other hand we do believe that the calculations described here achieve something new by 
dealing with the complete set of superficially divergent QCD vertices in one consistent approximation, 
and that they demonstrate a nontrivial, renormalization-related way of how the renormalization-group 
invariant mass scale establishes itself in the correlation functions of an asymptotically free theory. 
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Figure 1: DS equation for the r, 3 y Vertex. 
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Figure 2: Decomposition of the four-gluon T-matrix. 
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Figure 3: ’’Softened” exchange diagram (a) and diagrammatic representation of the compensating 
pole (b). 
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+ [terms not contributing to the 0-tli order] + ['2-loop terms] 


Figure 4: DS equation for the Bq part of the rsy-Vertex. 
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+ terms not contributing to 0-th order] 


Figure 5: DS equation for the ghost-gluon 3-point vertex in the ghost channel. 
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Figure 6: DS equation for the inverse fermion propagator. 
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+ terms not contributing to 0-th order] 



+ terms not contributing to 0-th order] 


Figure 7: Equivalent DS equations for the quark-antiquark-gluon vertex in (a) fermionic channel and 
(b) gluonic channel. 
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